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INTRODUCTION 


This  report  discusses  the  results  of  the  first  set  of  experiments  end  corre¬ 
sponding  theory  relating  to  the  evaluation  of  the  rise  tine  of  a  shock  front  in 
water.  This  report  also  discusses  the  planned  second  stage  of  experimentation 
and  corresponding  theory  wherein  outstanding  questions  are  to  be  addressed,  and 
where  a  preliminary  evaluation  of  the  shock  front  rise  in  the  explosive  liquid, 
nitromethane,  is  to  be  made. 

The  primary  reason  for  choosing  liquid  nitromethane  is  simple;  it  remains 
transparent  at  shock  pressures  below  those  necessary  for  shock  initiation  of 
detonation.  Furthermore,  experiments  have  demonstrated  a  shock-generated  voltage 
of  approximately  0.1  volt  (50  Q  termination)  in  liquid  nitromethane  at  50  kbar 
(ref  1).  If  that  0.1-volt  potential  drop  were  to  exist  entirely  across  a  small 
shock  front  thickness  (of  the  order  of  tens  of  molecular  diameters),  an  electri¬ 
cal  breakdown  mechanism  would  exist  for  bond  breaking  and  consequent  chemistry. 


THE  SHOCK  FRONT  THICKNESS  AT  5.8  kbar 

Two  recent  publications  (refs  2  and  3)  have  outlined  the  experimental  pro¬ 
gram  and  detailed  the  theory  leading  to  the  results  shown  in  figure  l.  This 
figure  indicates  that  a  combination  of  8^0.2  and  (An)f  _>  5.26  X  10-2  is  neces¬ 
sary  for  agreement  between  theory  and  experiment. 


The  calculations  illustrated  in  figure  1  were  carried  out  with  a  constant 
gradient. 


3n 

3x 


(An)f 

L~ 


constant,  n  “  index  of  refraction,  and  L 


shock  front  thickness 


/ 


shock  front  model.  8  =  (Lno/AQ),  and  (An)^  is  the  total  change  in  index  of  re¬ 


fraction  across  the  shock  front.  n„  is  the  unshocked  state  index,  and  A  is 

o  o 

optical  wavelength  in  vacuum. 


For  A  *  5145  A  and  nQ  «  1.33,  the  result  8  _<  0.2  gives  L  _<  623  A,  or  ap¬ 
proximately  215  molecular  lengths  for  an  oxygen-oxygen  distance  of  2.9  A.  (An)f 
=  5.26  X  10-2  corresponds  to  the  expected  index  of  refraction  change  at  5.8  kbar 
from  equation-of-state  (EOS)  data  (ref  3). 


With  the  constant  gradient  model,  the  reflectivity  It  can  be  put  in  the  form 

linO  \  2  ;  1/2 


R  *  sin* 


[ 1/2J’‘a  -  {  tj  [<*  +  1)2  {  1  -(; rrf)  *j 


co«9  -  (sin29  )ln 
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1  +  cos  9 
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DEGREES 


Figure  1 


.  Theory  and  experimental  data  for  shock  front  reflectivity 
at  5.8  kbar  in  water.  Curves  A,  B,  and  C  are  for 
1  optical  polarization  and  0  ■  (noL/Xfl)  ”0.1.  Curves 
D  through  G  are  for  II  optical  polarization  and 
(An)f  ”  5.26  x  10-2.  Constant  gradient  model 


calculations. 


(m  +l)2-2  sifl20 

1  -  tan2  0 


(m  +  l)2  -  sin20 


(2a) 


1  +  tan2  0 


(m  +  D: 


(m  +  l)2  -  sin20 


(2b) 


In  equation  1,  the  plus-and-mlnua  signs  refer  to  Incident  beam  optical  polariza¬ 
tion  perpendicular  and  parallel  to  the  plane  of  incidence,  respectively.  0  Is 
the  angle  of  incidence  in  the  unahocked  medium,  and  m  =  (An)/n  .  The  argument  of 
the  cosine  term  in  equation  1  represents  the  effect  of  the  opt?cal  phase  differ¬ 
ence  (at  the  reflectivity  measuring  photo-diode)  due  to  reflection  from  different 
layers  within  the  shock  front. 


While  the  constant  gradient  model  is  simple,  both  conceptually  and  algebrai¬ 


cally,  It  is  not  realistic;  the  foot  and  plateau  discontinuities  In 


3n 


3x 


(equlva— 


dp 

lent  to  discontinuities  in  P  “  mass  density)  would  not  be  expected  to  occur 
in  the  presence  of  viscosity.  A  more  realistic  (but  not  necessarily  correct) 
shock  front  model  is  the  hyperbolic  tangent  model  characteristic  of  weak  shock 
propagation  in  gases  (ref  4).  Indeed,  Cowan  and  Hornig  (ref  5)  have  employed  the 
hyperbolic  tangent  model  in  evaluating  shock  front  reflectivity  measurements  in 
gases. 


Writing  the  hyperbolic  tangent  model  as 


(An). 


1  +  e 


-4x7b 
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gives  for  the  reflectivity 
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The  difference  between  equations  1  and  4  is  entirely  due  to  phasing  effects.  The 
presence  of  the  m'  and  (ra^  -  ra"*)  terms  in  denominator  of  the  cosine  argument  of 
equation  4  is  representative  of  index-of-ref ractlon  changes  spatially  extending 
from  to  +“>. 

Although  algebraically  complicated,  equation  1  is  easier  to  Integrate  numer¬ 
ically  (by  partitions)  than  is  equation  4;  singularities  occur  in  the  inner  inte¬ 
gral  for  the  outer  Integral  limits.  If  the  mentioned  singularities  are  treated 
via  the  cutoff  procedure, 


*  10  ^  if  ra'<  m^  *  10  ^ ,  and 


C>) 


(mf  -  «')  =  af  x  10  ^  if  (m^  -  m')  <  x  10  \ 


(6) 


then  the  perpendicular  polarization  results  illustrated  in  figure  2  are  obtained 
via  numerical  integration  techniques.  The  figure  2  results  were  obtained  by 
using  (An)j  “  5.26  X  10-2  and  the  "DCADRE"  integration  program1  recursively. 
DCADRE  uses  a  Romberg  extrapolation  routine. 

Figure  2  shows  that  the  hyperbolic  tangent  results  are  essentially  the  same 
as  those  obtained  with  the  constant  gradient  model.  The  main  differences  between 
the  two  sets  of  results  are: 

1 • In  the  hyperbolic  tangent  case,  varying  0  •  Ln  /X  from  0.2  to  10-3 
no  longer  shifts  the  data  across  the  range  of  experimental ^at'a  (see  the  parallel 
polarization  data  and  calculations  of  fig.  1). 

2.  The  maximum  reflectivity  value  (that  obtained  for  0  m  10-3)  just 
coincides  with  the  lower  edge  of  the  data  band. 
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REFLECTIVITY 


HYPERBOLIC  TANGENT  CALCULATIONS 


OPTICAL  ANGLE  OF  INCIDENCE  (DEGREES) 

Figure  2.  Hyperbolic  tangent  shock  front  calculations  for  perpendicular 

optical  polarization  at  5.8  kbar.  Curves  A,  B,  and  C  are  constant 

gradient  model  results  (as  In  fig.  1).  All  hyperbolic  tangent 

-2 

results  are  for  (An).  “  5.26  x  10 


The  Implications  of  the  hyperbolic  tangent  model  are  essentially  the  same  as 
previously  reported  (ref  3)  for  the  constant  gradient  model:  first,  n  shock 
front  thickness  characterized  by  0  <  0.2  (equivalently  <215  molecular  lengths) 
and  second,  the  possible  necessity  of  a  non-equation-of-state  (non-EOS)  contribu¬ 
tion  to  the  final  Index-of-ref ractlon  change*  Such  a  non-EOS  contribution,  If 
positive,  would  Increase  the  calculated  reflectivity  for  a  given  value  of  0. 

For  a  given  value  of  0  and  (An)^,the  hyperbolic  tangent  reflectivity  will  be 
less  than  the  corresponding  constant  gradient  model  reflectivity  (because  An 
limbs  which  extend  spatially  from  -*  to  +“  Introduce  destructive  optical  Inter¬ 
ference).  As  0  Increases,  the  hyperbolic  tangent  destructive  Interference  In¬ 
creases  faster  than  Is  the  case  with  the  constant  gradient  model.  It  Is  precise¬ 
ly  this  optical  Interference  effect  which  prevents  the  hyperbolic  tangent  reflec¬ 
tivity  results  from  sweeping  across  the  experimental  data  band  as  0  Is  In¬ 
creased. 

Thus,  If  anything,  the  hyperbolic  tangent  model  results  reinforce  the  possi¬ 
ble  necessity  for  a  non-EOS  contribution  to  (An).  It  must  be  admitted,  however, 
that  the  necessity  for  introducing  a  non  EOS  contribution  to  (An)  is  not  glaring, 
and  a  definitive  conclusion  with  respect  to  the  (An)  origin  awaits  better  experi¬ 
mental  accuracy. 

It  is  difficult  to  measure  the  dimensions  of  an  object  with  a  comparison 
scale  (ruler)  which  Is  either  much  larger  or  much  smaller  than  the  object  to  be 
measured.  It  was  hoped  that  If  L«  A  was  Indeed  true,  the  spatially  extended 
(An)  limbs  associated  with  the  hyperbolic  tangent  model  would  effectively  yield 
a  ruler  and  I.  value  of  comparable  size.  Such  has  not  proved  to  be  the  case. 


SHOCK  FRONT  THICKNESS  PHYSICS 


The  question  of  shock  front  thickness  has  Intrigued  shock  wave  physicists 
for  many  years.  While  the  thickness  in  a  gas  is  well  understood  in  terms  of 
being  a  few  collision  mean  free  paths  (refs  5  and  6),  the  thickness  In  a  con¬ 
densed  medium  has  been  a  question  of  some  controversy  (e.g.,  the  mean  free  path 
of  what  excitation?). 

By  utilizing  computer  molecular  dynamics  (CMD),  techniques  D.  Tsai  et  al 
(refs  7  and  8)  have,  for  many  years,  demonstrated  results  which  are  consistent 
with  a  rise  time  occurring  over  a  few  lattice  parameters.  While  Tsai's 
results  were  Initially  for  solids,  other  workers  (refs  9  and  10)  had  reached 
similar  conclusions  for  monatomic  liquids.  Indeed  Tsai  et  al  (ref  11)  speculate 
that  5  to  10  lattice  parameters  (rise  time)  is  a  universal  effect  associated  with 
an  ’"Inertial  effect'  in  a  discrete  system." 

If  one  thinks  of  a  liquid  as  being  a  very  dense  gas,  and  as  such  not  con¬ 
taining  any  excitations  characteristic  of  long  range  order,  then  the  water  mole¬ 
cules  themselves  must  represent  the  only  entities  one  would  wish  to  consider  from 
a  mean-free  path  point  of  view.  Five  to  ten  lattice  parameters  (oxygen-oxygen 


n  i 


i 


6 


isstii 
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distances)  correspond  to  8  <  10  ,  a  number  not  Inconsistent  with  the  results 

demonstrated  in  figures  1  and  2. 

It  should  be  noted,  however,  that  to  the  extent  that  there  is  hydrogen  bond¬ 
ing  in  liquid  water  (refs  12  and  13)  there  is  the  possibility  that  such  hydrogen 
bonding  will  result  in  significant  long  range  order  and  that  the  long  range  order 
could  control  the  shock  front  rise  time.  At  this  point,  the  experimental-theo¬ 
retical  program  described  in  this  report  is  unable  to  decide  the  issue. 


IMPLICATIONS 


Normally  this  information  would  appear  at  or  near  the  end  of  a  report.  In 
this  case,  however,  the  implications  serve  at  least  as  a  partial  motivation  for 
the  following  two  sections  covering  the  next  set  of  theory  and  experiment. 

A  shock  front  thickness  of  less  than  30  A  would  have  both  basic  and  applied  i 
implications.  The  applied  aspect  pertains  to  the  possibility  of  a  new  concept  in  1 
the  shock  initiation  of  explosives,  while  the  basic  aspect  concerns  the  posslbil-  j 
ity  of  a  better  understanding  of  the  detailed  molecular  motions  occurring  within  j 
a  shock  front  in  a  liquid  (e.g.,  by  validating  the  CMD  approach).  j 

« 

As  previously  mentioned,  one  observes  a  potential  0.1-volt  drop  when  liquid  1 
nltromethane  is  shocked  to  near,  but  less  than,  50  kbars.  If  the  shock  front  j 
thickness  in  the  liquid  nltromethane  is  assumed  to  be  30  A,  and  if  the  potential  j 
drop  is  assumed  to  occur  only  across  that  shock  front  thickness,2  the  result  is  \ 
an  electric  field,  EgF,  within  the  shock  front  given  by 


ESF 


0.1  volts 
3  X  10~7  cm 


3  x  10'5  22!£b 

cm 


(7) 


The  value  of  Egp  given  by  equation  7  clearly  approaches  the  magnitude  necessary 
for  electrical  breakdown. 

If  electrical  breakdown  were  to  be  recognized  as  synonymous  with  bond  break¬ 
ing,  and  bond  breaking  synonymous  with  the  criteria  for  the  onset  of  chemistry, 
then  one  would  have  a  conceptually  simple  model  with  which  to  explain  shock  ini¬ 
tiation  of  explosives  —  a  shock  front  thickness  which  decreases  as  pressure 
increases  until  electrical  breakdown  and  chemistry  are  achieved. 


2While  it  is  well  known  that  many  liquids  have  an  electrochemical  ’’shock 
polarization”  contribution  (refs  14  and  15),  perhaps  a  sufficient  fraction  of 
the  0.1  volt  will  appear  across  the  shock  front  (due  to  the  orientation  of 
molecular  dipoles)  to  make  equation  7  essentially  correct. 


7 


The  two  Ingredients  associated  with  the  electrical  breakdown  concept  are  the 
shock-induced  "polarization"  voltage  and  the  front  thickness  as  a  function  of 
shock  amplitude  (pressure).  These  two  Ingredients  serve  as  the  basis  for  the 
next  set  of  theory  and  experiments.  The  front  thickness  as  a  function  of  pres¬ 
sure  will  be  pursued,  as  before,  from  a  continuum3 *  optical-mechanical  point  of 
view.  The  polarization  voltage,  however,  In  sofaras  It  Involves  the  partial 
alignment  of  molecular  electrical  dipoles,  is  a  molecular  level  calculation. 

Our  recent  publications  (refs  2  and  3)  attempted  to  show  a  relationship 
between  a  5.8  kbar  shock  front  thickness  upper  limit  of  approximately  215  oxygen- 
oxygen  distances  and  Soviet-reported  anomalously  high  viscosity  values  of  approx¬ 
imately  101*  P  in  the  100-kbar  region.  Unfortunately,  the  hyperbolic  tangent 
results  illustrated  In  figure  2  maintain  the  previously  reported  shock  front 
thickness  upper  limit  (refs  2  and  3)  and  thus  do  not  resolve  the  viscosity  ques¬ 
tion.  It  Is  worth  noting,  however,  that  Hamann1*  (ref  17)  casts  doubt  on  the 
Interpretation  of  some  of  those  Soviet  shock  viscosity  experiments  and  also 
doubts  that  some  of  the  Soviet  experiments  have  measured  shock  front  (as  distinct 
from  behind  the  shock)  viscosity  as  we  had  previously  maintained  (ref  3). 

One  possible  explanation  for  the  Soviet  viscosity  observations  could  be  the 
breaking5  of  hydrogen  bonds  in  the  100-kbar  region.  As  hydrogen  bonds  result  In 
a  restoring  force  which  acts  to  counter  a  net  alignment  of  molecular  dipoles  (see 
the  section  on  Shock  Polarization),  bond  breaking  should  increase  observed  shock 
polarization  voltage  amplitudes  for  water.  This  explanation  Is  consistent  with 
the  concept  by  Graham  (ref  18)  wherein  bond  breaking  Is  claimed  to  be  responsible 
for  the  observed  shock  polarizability  of  polymeric  solids.  Thus,  the  planned 
shock  polarization  experiments  and  their  theoretical  interpretation  should  be 
vary  pertinent  to  the  question  of  Soviet-claimed  viscosity  values. 


REFLECTIVITY  VERSUS  PRESSURE 


The  upcoming  ref lectlvlty-versus-pressure  experiments  will  be  analyzed,  as 
was  done  previously,  by  using  the  constant  gradient  and  hyperbolic  tangent  models 
for  data  at  each  pressure  point.  Additionally,  there  Is  Information  to  be  gained 
from  any  change  in  the  reflectivity  as  a  function  of  pressure. 


3Water  would  be  a  very  interesting  material  to  treat  via  CMD  shock  calcu¬ 
lations.  CMD  (static)  calculations  have  been  carried  out  (ref  16)  yielding 
an  "effective"  intermolecular  potential,  and  that  intermolecular  potential  is 
available  for  the  CMD  shock  calculations. 

‘‘Private  conversation. 

5Slnce  bond  breaking  requires  energy,  such  an  energy  loss  could  be  Interpreted  as 
a  viscous  phenomenon. 
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By  taking  the  derivative  of  equation  1  with  respect  to  pressure,  it 
shown  that 
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of  mj 


At  5.8  kbar  one  has,  from  EOS  considerations,  -  0.0405.  For  that 
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can  be 
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value 
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The  term  in  equation  8a  containing  the  derivatives  of  lnmj  and  InL  for 
(n  L/X  )  =  10  1  is  smaller6  by  a  factor  of  approximately  0.2  than  the  first 


Thus,  there  should  be  sufficient  accura- 


terra,  which  is  evaluated  at  (6^,  mf). 

3R 

cy  in  the  difference  between  -5—  and  the  first  term  on  the  right  of  equation  8a 

3p  f  31nmf  31nL^ 

to  determine  at  least  the  sign  of  {  — -  Sp~  J  *  hydrogen  bond  break¬ 

ing  were  to  lead  to  a  large  increase  in  L  with  increasing  shock  pressure,  that 
increase  in  L  would  become  obvious  from  sign  considerations  alone. 


SHOCK  POLARIZATION  CONSIDERATIONS 


In  this  calculation  the  net  orientation  of  a  hydrogen-bonded  water,  based  on 
the  "rigid  four-point-charge"  model  molecule  concept  (ref  9),  is  crudely  calcu¬ 
lated.  In  the  static  CMD  simulation  calculations  (ref  19),  each  four-point-char¬ 
ge  molecule  interacts  via  an  elecrostatlc  potential  function.  Because  the  calcu¬ 
lation  presented  below  considers  a  single  such  molecule  in  an  average  mechanical 
force  field,  the  inertial  effect  of  hydrogen  bonding  is  crudely  represented  by 
simply  placing  a  hydrogen  atom  mass  at  each  of  the  four  charge  points.  See 
figure  3. 

By  employing  Euler  angle  transformations  (ref  20)  and  a  net  z-directed  force 
(due  to  the  shock)  on  each  hydrogen  atom,  it  is  possible  to  derive  an  equation  of 
motion  for  0  .  Using  the  Euler  angle  geometry  shown  In  figure  4  allows  0  to  be 
written  as 


na0H9 


21  ma  f 
o  OH  n 


JTT 


sinQ 


-  -  ( 
2  ^ 


OH 


OH 


(9) 


where 


OH 


”  oxygen-hydrogen  distance  ~  1A 


*HB 


«  hydrogen  bonding  distance  ~  2A 


6Because  the  argument  of  sine  term  is  small, 
approximately 


277 

m 

f 


n  L 
o 

X 


V 

J 


the  total  factor  to  be  evaluated  is 
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Figure  3.  The  tetrahedral  geometry  of  the  rigid  four-point-charge  model. 


Figure  4 


Pq  is  the  net  molecular  polarization,  and  the  are  the 
coordinate  vectors  of  the  charges  (and  associated  hydrogen 
masses). 


z 


The  Euler  angle  coordinates.  The  rigid  coordinate  system 
attached  to  the  rigid  molecule  of  figure  3  becomes  the 
(x't  y",  z")  system  and  is  allowed  to  rotate  within  the  fixed 
system  (x,  y,  s).  (0,  4,  f)  are  Euler  angles. 


hydrogen  atom  mass 


^  * 


m 

*o  ■  N  -  |M  “  °-2A 

Ig  »  moment  of  Inertia  about  oxygen  atom  center  with  respect 

4  „ 

to  line  of  nodes  axis  ■  y«  (1.64)  (a0jj)  ■ 

fn  »  average  force  on  each  hydrogen  atom  *  a  • 

s 

k  -  elastic  restoring  force  coefficient 

5$1  =*  0  -  0q 


Equation  9  is  the  three  dimensional  version  of  a  calculation  first  presented 
by  Horle  (ref  21).  In  f  ,  a  is  the  fraction  of  the  force  per  water  molecule 

acting  on  each  hydrogen  mass,  a  -  1/18,  a  is  the  shock  stress  amplitude,  N  the 

number  density  of  hydrogen  molecules,  and  L  is  the  shock  front  thickness.  It 

s 

Ip  assumed  that  fn  acts  only  within  the  shock  front.  The  term  arises  from 

the  stretching  of  the  hydrogen  bond  acting  as  a  restoring  force  to  the  angular 
motion. 

The  restoring  force  derivation  proceeds  with  the  aid  of  figure  5  and  the 
assumption  of  small  angle  changes,  as  follows. 


®HB 


f  («2)21 

•m  C086®2  "  I1 "  “T  ) 


(10) 


(«e2)2 

"  aHB  ”  ^IB  ~2~  ‘ 


(11) 


(«e2)2 

/,  F  -  -  k(a^B  -  aHR)  -  "  k  am  -y-  . 


(12) 


1?. 


(13) 


(6&2>3 

6F  “  F  sin  602  “  "  k  — j-  • 


-us) 


■ro,  <«M3- 


Setting  60i  »  0  -  0q,  where  9o  denotes  the  unshocked  0  coordinate  (thermal  mo¬ 
tions  being  neglected)  of  a  particular  molecule,  is  consistent  with  the  crude 
calculation  presented  here. 

Equation  9  allows  two  interesting  limiting  cases:  setting  k  *  0  describes 

the  free  (bond  broken)  molecule,  or  0  “  0  describes  the  hydrogen  bond  controlled 
restoring  force  limiting  angle. 

From  equation  9  with  0*0, 


164  NL  )  Sln8 

6 '  °°)3  ’  ' 


o  -  0  + 

O 


1 / 


.  l/3« 

sin  0. 


With  the  coefficient  of  the  sin  0  term  assumed  small 


0  -  0  + 


( 7t.  -  «Y 

V  ML  jl*  / 


,  1/3- 
sin  0o, 
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(18) 


<CO80> 


-  1/2 


7 A  a  al1'3 

-L- J  J 

“VbiJ 


/'sin 
1  o 


7/3 


e  de  , 

o  o’ 


(19) 


where  <cos0>  denotes  the  average  value  obtained  by  Integrating  cosO  over  all 
possible  0q  positions  (on  a  unit  sphere). 

Carrying  out  the  Integration  of  equation  19  by  simple  numerical  techniques 

gives 


<COS0> 


(0.74) 


‘7 lQ  a  o']173 

-  •‘"VohI  ’ 


(20) 


Employing  k  -  w  2M,  M  being  the  water  molecule  mass  (the  rigid  water  molecule 
being  hydrogen  bonded  to  the  surrounding  water  structure),  with  <*>  corresponding 
to  approximately  9  kcal  per  mole  (ref  22),  along  with  N  »  3.47  x  1022  cm-3, 

a  »  1010  cgs  (i.e.,  10  kbars),  and  L  -  6.2  x  10  6  cm  (corresponding  to  0  »  0.2 
In  f igs  1  and  2)  gives  8 


<cos0>  *  -  2.4  x  10"4. 


(21) 


The  net  shock  front  produced  polarization  (in  the  shock  direction)  Is  thus 
given  by 


IP  -  -  HP  <CO80>, 
o  ’ 


(22) 
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esu  cm  from  figure  3 


or  with  P 

o 


2.27  x  10"18 


T  -  1.9  x  10  cgs7  . 


(23) 


The  weak  shock  current  density  at  t  "  0  (when  the  shock  first  enters  the 
water)  is  given  (ref  23)  by 


Jio 


TO 

L  * 
o 


(24) 


with  U  being  the  shock  velocity,  L  the  water  sample  thickness,  and  the  subscript 
10  denoting  10  kbars.  For  U  =*  2  x^O5  cgs,  and  Lq  “  1  cm, 


j I  0  *  3.8  x  106  cgs. 


(25) 


Converting  to  MKS  units  via  J\q  “  jio  /(12rr  x  10s),  and  assuming  a  sample  cross 
section  of  10--4  mz,  finally  gives  for  the  weak  shock  polarization  current  in 
water 


III)  -  10 


Amperes. 


(26) 


Equation  26  predicts  a  10  kbar  shock  polarization  voltage  across  a  50S  term¬ 
ination  load  of  5  mv.  Evaluating  equations  20  through  24  at  100  kbar.  with  U 
taken  as  5  x  105  cgs  (ref  23)  and  all  other  parameters  unchanged,  yields  25  mv 
for  the  polarization  signal  across  a  500  termination.  That  25-mv  result  is  too 
small  by  one  order  of  magnitude  when  compared  with  the  assumed  100-kbar  observa¬ 
tions  of  Eichelberger  and  Hauver  (ref  24). 

There  are  a  number  of  possibilities  for  the  small  25-mv  result: 

1.  is  approximately  two  orders  of  magnitude  smaller  than  the  upper 
limit  value  of  6.2  x  10-6  cm  employed  above. 

2.  The  free  (i.e.  k  ”  0)  molecule  version  of  equation  9  is  making  at 
least  a  partial  contribution  to  IP  . 


'The  actual  molecular  dipole  moment  is  a  complicated  function  of  the  hydrogen 
bonding  in  liquid  water.  (See,  for  example,  ref  22.) 


3.  Little  of  this  section  corresponds  to  reality. 


Because  the  free  molecule  (i.e.,  k  ■  0)  version  of  equation  9  can  be  shown 
to  yield  <cos6>  results  in  orders  of  magnitude  larger  than  those  associated  with 
equation  20,  the  25-mv  calculated  value  can  be  seen  as  supporting  the  hypothesis 
that  significant  bond  breaking  occurs  at  presaures  approaching  100  kbars.  Such 
bond  breaking  would  add  a  pressure  dependent,  free-molecule  contribution  to  IP  , 
and  would  be  consistent  with  ideas  put  forward  by  Graham  (ref  18). 


STATUS  SIMMARY 


The  question  of  shock  front  thickness  in  water  remains  unanswered  except  for 
a  predicted  upper  limit,  at  3.8  kbar,  of  approximately  6  x  10-6  cm.  Similarly, 
the  anomalmously  large,  Soviet  shock-related,  viscosity  values  near  100  kbar  are 
neither  denied  nor  confirmed  by  the  theory  and  experiment  reported  herein. 

Extensive  hydrogen  bond  breaking  in  liquid  water  could  be  an  explanation  for 
some  guess  extrapolations  to  pressures  above  5.8  kbars;  however  the  work  carried 
out  to  date  (October  1981)  is  Insufficient  to  make  a  positive  judgment  on  this 
point. 

Hopefully,  future  theory  and  experiment  outlined  in  this  report  will 
decisively  settle  all  of  the  outstanding  questions. 
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